Chapter 5 Application of Bi - CGSTAB to discretized coupled PDEs
نویسنده
چکیده
A version of the Bi-CGSTAB method is applied for solving the linear systems that typically occur when applying the Newton method to a discretized set of coupled elliptic partial di erential equations in two dimensions. The Incomplete Line LU decomposition is generalized for the case of coupled equations and applied as a preconditioner. A suitable stopping criterion is developed for the Bi-CGSTAB method. Numerical experiments are performed for problems stemming from the area of semiconductor modelling and for an aquifer problem. For the latter problem a comparison is made with an existing multigrid algorithm in which Incomplete Line LU relaxation is used as smoothing procedure. Note: This chapter is a revision of Report NM-R9213, CWI, Amsterdam, 1992. 5.
منابع مشابه
Bicgstab(l) for Linear Equations Involving Unsymmetric Matrices with Complex Spectrum
For a number of linear systems of equations arising from realistic problems, using the Bi-CGSTAB algorithm of van der Vorst [17] to solve these equations is very attractive. Unfortunately, for a large class of equations, where, for instance, Bi-CG performs well, the convergence of BiCGSTAB stagnates. This was observed specifically in case of discretized advection dominated PDE’s. The stagnation...
متن کاملBicgstab(l) for Linear Equations Involving Unsymmetric Matrices with Complex
For a number of linear systems of equations arising from realistic problems, using the Bi-CGSTAB algorithm of van der Vorst [17] to solve these equations is very attractive. Unfortunately, for a large class of equations, where, for instance, Bi-CG performs well, the convergence of BiCGSTAB stagnates. This was observed specifically in case of discretized advection dominated PDE’s. The stagnation...
متن کاملMultigrid Preconditioners for Bi-cgstab for the Sparse-grid Solution of High-dimensional Anisotropic Diffusion Equation
Robust and efficient solution techniques are developed for high-dimensional parabolic partial differential equations (PDEs). Presented is a solver based on the Krylov subspace method Bi-CGSTAB preconditioned with d-multigrid. Developing the perfect multigrid method, as a stand-alone solver for a single problem discretized on a particular grid, often requires a lot of optimal tuning and expert i...
متن کاملAdaptive moving mesh computations for reaction-diffusion systems
In this paper we describe an adaptive moving mesh technique and its application to reaction-diffusion models from chemistry. The method is based on a coordinate transformation between physical and computational coordinates. The transformation can be viewed as a solution of adaptive mesh partial differential equations (PDEs) which is derived from the minimization of a mesh-energy integral. For a...
متن کاملNumerical Study of Iterative Methods for the Solution of the Dirichlet-Neumann Map for Linear Elliptic PDEs on Regular Polygon Domains
A generalized Dirichlet to Neumann map is one of the main aspects characterizing a recently introduced method for analyzing linear elliptic PDEs, through which it became possible to couple known and unknown components of the solution on the boundary of the domain without solving on its interior. For its numerical solution, a well conditioned quadratically convergent sine-Collocation method was ...
متن کامل